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ABSTRACT 


The phenomenon of dynamic fracture is of interest in a wide variety of contexts from 
high-speed manufacturing and defence applications to geophysics ^d seismology. In 
dynamic fracture problems, role of material inertia becomes significant The formulation 
and analytical solutions of dynamic fracture problems are very complex. Experimental 
and numerical methods like finite element method are used extensively for solving 
dynamic fracture problems. Existing methods for simulating crack phenomenon in 
dynamic fracture problems such as, moving mesh procedures or force release models, 
have certain drawbacks that they are either computationally very intensive or give 
oscillations in the solution. 

Present work is a modification of ‘Stiffness Release Model’ [19, 20] developed 
for simulating high-speed dynamic crack propagation. In the ‘Stiffiiess Release Model’ 
the mesh is kept stationary and a 1 -Dimensional elastic spring element is added at the 
crack tip node. The gradual propagation of the crack from one node to the next node is 
achieved by decreasing the stiffhess value from a very large value (ideally infinite) to 
zero as the crack propagates to the next node within the element. But this model gives 
high oscillation at the element boundaries. 

The emphasis of the present work is to improve the above ‘Stiffiiess Release 
Model’ by mass modification. Here mass of the element surrounding the crack tip is 
modified so that it varies as a function of non-dimensional crack length when crack 
extends through the element and the effect of other parameters with crack velocity is also 
studied which will enable in finding the propagation history under impact loading. 

The finite element model is first validated for static case and then with 
hypothetical problem for different crack speeds. Results give smooth variation of energy 
release rate curve for different crack velocities. 
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CHAPTER 1 


INTRODUCTION 


1.1 Introduction 

The phenomenon of dynamic ftacture is of interest in a wide variety of contexts from 
high-speed manufecturing and defence applications to geophysics and seismology. A 
considerable amount of research has been directed towards the solution of problems 
involving the interactions of stress waves with cracks and boundaries to improve 
understanding of the behaviour of material failure under dynamic loading. The stu(fy 
of crack interaction is also of great interest from the point of view of quantitative non- 
destructive evaluation of structural materials. 

The (fynamic fracture phenomenon can be characterised by various cfynamic 
states of crack tip. The dynamic states of crack tip are induced by impact loa ding 
£q)plied to the cracked soHds, or by fiast motions of the crack tip itself Dynamic loads 
may be created by moving vehicles, wind, gusts, shocks or blasts, unbalanced machines, 
wave impacts, seismic disturbances, etc. For solving these jxoblons, oftei the time 
variation is ne^ected and they are treated as static or quasi-static problems. Structural 
members, which are subjected to cychc loading such as vibrations, microscopic 
imperfections in the material give rise to microscopic cracks which in turn grow into 
cracks of significant size over passage of time. At this stage, further growth of cracks 
stably or unstably leading to the component feilure depends not only on the material 
characteristics but also on the nature of loading, i.e., static or dynamic. For many 
dynamic problems, it is impossible to find closed form solutions, especially if they 
involve the aspect of fracture. A number of analytical solutions to problems of (fynamic 
j&acture shpdding important h^ on the basic phermienon have been rqwrted in the 
past two decartes. These solutions are however limited to simple cases of loading and of 
unbounded plane, bodies. Usually, the interactions of the stress waves emanating finm 
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the crack tip and those reflected fix>in the boundaries of a finit e body make the analytical 
solutions invahd. Thus, it often becomes mandatoiy to analyse dynamic crack 
propagation in finite sohds using numerical methods. 


1.2 DyBamic Fracture Mechanics 

Dynamic fracture mechanics is the subfield of firacture mechanics concerned with 
fracture phenomena for which the role of mate rial inertia becomes si gnifi cant One 
method of classifying dynamic fracture problems is as follows (Nishoka and Atulri, 
1986, [1]); 

(1) Solids containing stationary cracks subjected to (fynamic loading. 

(2) Solids containing dynamically propagating cracks under quasi-static 
loading. 

(3) Solids containing dynamically propagating cracks under dynamic loading. 

A knowledge of time-dependent asymptotic stress and displacement fields 

near the crack tip, and their "strength" as quantified, for instance, by the "stress 
intensity factor (SEF) in linearly elasto-dynamics, is basic in xmderstanding the 
process and nature of dynamic firacture in solids. A direct apjnoach to laboratory 
evaluation of dynamic firacture toughness and crack arrest toughness, is based on the 
laboratory measurements of the instruments dynamic stress fields close to the 
propagating crack tip using photo-elastic method. Such direct measurements, 
however, are generally difficult To overcome these difficulties, and simplify the 
measurements necessary, hybrid experimental-numerical methods are often preferred 
This is one reason that the advancement of the state of science of dynamic fincture 
mechanics relies heavily on simultaneous advances in computational methods. 

Dynamic fincture mechanics concerns (i) the onset of crack growth under 
dynamic loading (ii) (fynamic crack popagation in stressed solids. The conventional 
linear elasto-plastic or quasi-static fincture mechanics applies only up to the end of 
stable growth under quasi-static loading and assumes that the onset of unstable crack 
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propagation renders the structure useless. Yet, the prevention of crack growth 
initiation itself may be too conservative or too costly an objective for some structural 
designs, and also catastrophic failures caused by unstable crack growth are obviously 
intolerable. In these cases, the assurance of crack arrest, as a second lin e of defence, 
is essential. This concept has been investigated for design methodologies assuring 
the integrity of nuclear pressure vessels under ther mal shock conditions, LNG ship 
hulls, and gas transmission pipelines. 

The fundamental flame work of the subject of dynamic fiacture mechanics, 
relies primarily on solutions for dynamic behaviour of solid containing cracks these 
solutions are characterized by, (i) stress wave integration (ii) the need to account for 
kinetic energy in the global energy balance of the flacturing body, and (iii) inertia 
effects of the material. The inertia effect may arise due to two reasons (i) rapidly 
applied load on a cracked solid and (ii) rapid crack propagation. In the first case the 
influence of the load is transferred to the crack by stress waves throu^ the material. 
In the second case, the material particles on the two crack feces displace with respect 
to each other, as the crack advances. In the first case, inertia effect will be significant 
if the characteristic time (i.e. maximum load by the rate of load increase) is larger 
than time required for stress wave to travel at the characteristic wave speed of the 
material over a representative length of the body (crack length or distance fiom crack 
edge to the loaded boundary) (Freund, 1990, [2]). In the second case, the inertia effect 
should be accounted for whenever the crack velocity is a significant fiaction of the 
characteristic wave speed of the material. 

In the past two decades or so, a number of analytical solutions, which provide 
a useful understanding of dynamic crack behaviour, have been obtained. These 
analytical solutions are however, limited to cases of simple loading and unbounded 
plane bodies. Moreover, the stress wave introductions, which play an important role 
in dynamic fiacture mechanics, usually render the analytical solutions intractable. 


3 



Therefore, the use of numerical methods is often indispensable for the analysis of 
cracks in finite solids. 

The fimte element method was more suitable for the analysis of stationary 
cracks under dynamic loading, due to the fact that the relevant singularities can be 
modelled in the crack tip elements. In large measure, this is the result of the 
development of finite element procedures to model propagating singularities and of 
path-independent integrals that characterize the strength of the fields near prop agating 
crack tips. 

13 Energy Concepts in dynamic fracture mechanics 

Analytical methods based on the work done by applied loads and the changes in the 
energy of a system that accompany a real or virtual crack advance have been of 
central importance in the development of dynamic fincture mechanics. Griffith (1920, 
[3]) recognized the importance of the variation of the energy measures during crack 
growthL Gjnsequently, for an elastic body con taining a crack, the negative of the rate 
of change of total potential energy with respect to crack dimension is called the 
‘Energy Release Rate’, G and is a function of crack size in general. G is the amount 
of energy, per unit length along the crack edge that is suppUed by the elastic oiergy in 
the body and by the loading system in creating the new fiacture surftice. 

13.1 Dynamic Energy Release Rate 

The cfynamic energy release rate is defined as the rate of mechanical energy flow out 
of the body and into the crack tip per unit crack advance. The energy flux, or energy 
flow per unit time, must be divided by v, the crack movement per unit time, in order 
to obtain the energy release from the body per unit crack advance. For two- 
dimensional fields, G is the mechanical energy released pa: unit thickness of the body 
in the direction perpendicular to the plane of deformatiotL 
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1.4 Rapid (Fast) Crack Propagation 

Analyses of fest crack growth in structural mechanics have generally been based on 
an approach where a material and crack speed dependent value of a characterizing 
parameter, such as the energy release rate or the stress intensity factor, is used in 
conjunction with a crack tip equation of motion, Freund(1990,[2]). 

1.4.1 Fast crack growth at constant speed 

The assumption of constant speed results in significant simplification of the boundary 
value problem to be solved for the governing partial differential equations. The 
restriction to the constant speed crack growth permits a reduction of the number of 
independent variables by one and the analysis is simplified considerably. But some of 
the features of the crack growth process that are physically significant are overlooked 
in this approach. 

1.4.2 Crack growth at nonimiform speed 

The assumption of constant speed growth is imposed in order to render the 
mathematical models tractable. The study of a problem involving crack growth at 
nonuniform speed proceeds in two steps 

1) The boundary value problem is considered for arbitrary motion of the crack 
tip, with a view toward obtaining a full description of the mechanical fields 
near the crack edge during growth. 

2) An addition physical postulate in the form of a crack propagation criterion is 
imposed on the mechanical fields in order to determine an equation of motion 
for crack tip. The analytical approaches developed for stationary cracks and 
cracks growing at a constant speed can’t extend to the case of nonuniform 
crack growth. 
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1.5 Fracture Model for dynamic crack propagation 

Over the last few years, physicists interested in nonlinear phenomena have returned to 
problems of crack propagation, and to the study of model problems in which cracks 
develop from the breaking of bonds. Also, it is now possible to model the entire 
development of a crack, in a finite-element model of a continuum, even with inelastic 
effects (that is, plasticity) allowed for. Damage mechanics is well developed at the 
phenomenological level, and can also be built into codes that predict fracture. The 
admission of any model for microscopic processes introduces a length scale, and size 
effects may emerge. These are of pjractical importance as well as theoretical interest 
Mathematical theory based on variatioiial methods is developing, at least for cracks in 
elastic media. However, many problems remain: there is no adequate theory for the 
stability of a propwgatir^ crack, even in the framework of linear elasticity, at the 
present time. Interactions between macro- and micro-cracks require better modelling. 
The developjment of (approximate) fractal structure has so far not been modelled 
mechanistically. Some hint of the influence of underlyit^ lattice structure has been 
demonstrated (it may ap)p»arently prevent steady-state propagation at low sp)eed) but 
this is not fully wo±ed out 


1.6 Literature Survey 

Though considerable amount of work has been done to study the fracture phenomena 
through numerical methods, only rectangular double cantilever beam (DCB) 
sp)ecimens were given sp)ecial interest. Owen and Shantaram (1977, [4]) started the 
case of finite element method to stucfy dynamic crack growth. They studied DCB 
sp)ecimen and pipeline problem under transient loading. The crack was advanced from 
one node to the next node and when stress at the gauss px)int nearest to the crack tip 
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exceeded a certain value, damping coefficient was made zero at the released node. 
The crack propagation history was simulated for the given loading conditions. 

Nishioka and Atluri (1982b, [5]) investigated the crack propagation and arrest 
in high strength steel DCB specimen using moving singular dynamic finite element 
procedme. An edge crack in a rectangular DCB specimen was propagated by 
inserting a wedge and the results were compared with e)q)erimental data In another 
work, Nishioka and Atluri (1982a, [6]) presented the results of generation and 
prediction studies of dynamic crack propagation in plane stress and plane strain cases. 
The studies were conducted by using FEM, taking into account stress singularity near 
the crack tip. Variation of dynamic stress intensity fector with time and variation of 
dynamic fincture tou^huess with velocity were studied and compared with available 
experimental results. 

Nishioka and Atluri (1986, [1]) gave elaborate information on the analysis of 
dynamic firacture using FEM. The most common way to deal with the crack tip region 
was to simulate crack growth through gradual release of elemental nodal forces or 
embedding a moving element in ^^1lich interpolation fimctions are determined by the 
continuum near tip fields at the crack tip in the mesh, and J-Integral consideration. 
Following special description the differential equation in time for the node point 
variables must be integrated. Because the (fynamic fields associated with rapid crack 
growth are rich in high ficquency content, smaU time steps are required for accuracy. 
Because of natural restriction to small time steps, the author s report that many times 
combination of conditionally stable explicit time integration scheme and diagonalised 
mass matrix are found to be accurate. 

Chiang (1990, [7]) presented a numerical procedure based on eigen function 
to determine the dynamic stress intensity factor for a crack moving at steady state 
under antiplane strain condition. An edge crack problem and a radial crack problem 
were solved using traction and displacement boundary conditions separately. It was 
shown that the dynamic effect was relatively insignificant for low crack propagation 
speeds provided that specimen size was fairly large. 
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Thesken and Gudmundson (1991, [8]) worked with elasto-dynamic moving 
element formulation incorporating a variable order singular element to enhance the 
local crack tip description. Kennedy and Kim (1993, [9]) incorporated micropolar 
elasticity theory into a plane shain fimte element formulation to analyse the cfynamic 
response of the crack. Materials with strong micropolar properties were found to have 
sigmficantly lower dynamic energy release rate than their classical material 
counterparts. 

Wang and Williams (1994, [10]) investigated high-speed crack growth in a 
thin double cantilever beam specimen using FEM. The cantilever end was loaded 
under an initial step displacement of 1mm and then pulled with a constant velocity the 
crack was propagated at assumed speed by gradual node release technique. Large 
dynamic effects were observed because of wave reflections within fini te specimen 
size. The reflection and interaction of the waves from free boundary were avoided by 
limiting the duration of study. Chandra and Krauthammer (1995, [11]) investigated 
effects of rapidly applied load on the J-integral and stress intensity factor (K) of a 
cracked solid. Variation of J has been studied using energy balance apparoach, 
wbereas that of K has been dealt with elastodynamic considerations together with 
several simplifying assumptions and approximations. Lee, Hawong and Choi (1996, 
[12]) studied propagating crack problems of orthotropic material under the dynamic 
plane mode. Dynamic stress components and dynamic displacement components 
around the crack tip of an orthotropic material under the dynamic load and the steady 
state in crack propagation were derived. When the crack propagation speed 
approaches zero, dynamic stress components and dynamic displacement components 
are identical to those of static state. The stress component values of the crack tip are 
greater when the fiber direction coincides with the direction of the stress component 
than when the fiber direction is normal to the stress component 

Lin and Smith (1997, [13]) presented a method wherein crack growto is 
predicted in a step by step basis from Paris law using stress intensity factor calculated 
by using finite element method. The crack front is defined by a cubic spline curve 
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from a set of nodes. Both the one quarter node crack opening and 3-D J-Integral 
method are used to calculate stress intensity factor. Automatic remeshing of finite 
element model to a new position that defines the new crack front enables the crack 
propagation to be followed Beissel, Johnson and Popelar (1998, [14]) presented an 
algorithm, which allows crack propagation in any direction but doesn’t require 
remeshing or the definition of new contact surfeces. This is achieved by tracking the 
path of the crack tip and failing the element crossed by the path such that they can no 
longer sustain tensile volumetric stress. The edges of these feiled elements simulate 
crack faces that can sustain only compressive normal tractiotL 

Zhnang and Guo (1999, [15]) addressed recent developments in the area of 
dynamic fracture mechanics and the apphcations of analysis methods to the rapid 
crack propagation for gas pipelines. Criteria for crack initiation, propagation and 
arrest were discussed Christina and Christer (2001, [16]) presented a method for 
obtaining the complex stress intensity factor for an interface crack in a bimaterial 
using minimum number of computations. A crack closure integral method has been 
used Falk, Needleman and Rice (2001, [17]) performed simulation on a square block 
in plain strain with an initial edge crack loaded at a constant rate of strain. 

Brener and Spatschek (2002, [18]) presented a continuum theory which 
describes the fast growth of a crack by surface diffusion. It predicts the saturation of 
steady state crack velocity appreciably below the Raylei^ speed and tip blunting. 


9 



1.4 Present Work 

Present work is a modification of ‘Stiffiiess Release Model’ [19, 20] developed for 
simulating high-speed dynamic crack propagation. In the ‘Stiffiiess Release Model’ 
developed by Kishore, Reddy and Deshmukh [19,20], the mesh is kept stationary and 
a 1 -Dimensional elastic spring element is added at the crack tip node. The gradual 
propagation of the crack from one node to the next node is achieved by decreasing the 
stif&iess value ft^om a very large value (ideally infinite) to zero as the crack 
propagates to the next node within the element. But this model gives high oscillation 
at the element boundaries. 

The emphasis of the present work is to improve the above stiffiiess Release 
Model by mass modification. Here mass of the element surrounding the crack tip is 
modified so that it varies as a fimction of non-dimensional crack length \\hen crack 
extends through the element and the effect of other parameters with crack velocity is 
also studied which will enable in finding the propagation history under impact 
loading. 

Chapter 2 describes the basics of FEM formulation and the existing 
propagation model. 

Chapter 3 presents the crack propagation element model. 

Chapter 4 describes the results and discussion 

Chapter 5 presents the conclusion of the present work and scope of future 

work. 
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CHAPTER 2 


FEM FORMULATION AND CRACK MODELEVG 

2.1 Basics of Dynamic Analysis of Plates with Crack 

Plates with cracks subjected to dynamic load is solved using the theory of flexural motion of 
plates developed by Mindlin (1951, [21]) in which the three physical boundary conditions, 
viz. applied bending moments, twisting moment and transverse shear can be satisfied 
independently at the crack edge, whereas, the approximate Kirchhoff method can satisfy only 
two boundary conditions on the crack face leading to unrealistic results. 

2.1.1 Mindlin Plate Theory 

Consider the plate shown in Figure 2.1. The X and Y-axes are defined to be in the mid plane 
of the plate and z-axis is normal to the plane of the plate. The corresponding displacement 
components are u, v and w, respectively. Making the following assumptions; 

1) The strains and stresses in the Z direction are negligible 

2) The normal to the middle plane of the plate remains straight during deformation 

static equilibrium equations for the plate bending problem can be obtained (Zienkiewicz, 
1991, [22]; Timoshenko, (1959, [23]) in terms of mid plane rotations (0x, 0y), lateral 
displacement (w) of the mid plane and shear force (S) at the mid plane. The representations 
of variables used in plate approximation are shown in Figure 2.2. 

For dynamic analysis inertia effects have to be considered Using D’Alemberts 
principle of motion, forces due to mass and moment of inertia, are taken into account m the 
equilibrium equations. With reference to the convention shown in Figure 2.2 the final 
equations are modified as follows ( Min dl in , 1951, [21]). 
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Figure 2.1 Cracked plate subjected to lateral dynamic load 
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p is the mass density of the plate. 
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M ^ and are moments on faces perpendicular to X and Y axes due to anda„ 

stresses respectively, 

is the twisting moment due to shear stress , 

0^ and 6y are the rotation of faces perpendicular to the X and Y-axes respectively, 
and Sy are shear forces on faces perpendicular to the X and Y-axes respectively, 
and Sy are strains along X and Y-axes respectively, and 
andYyj are shear strains in X-Z and Y-Z planes respectively. 


2.2 Finite Element Formulation 
2.2.1 FEM Equation 

The problem is formulated in the three sets of variables 0, S and w, which are expressed in 
terms of nodal parameters using appropriate shape functions and parameters as 


0 = 



= Ne0 


S = 



( 2 . 6 ) 


w = N^w 

where Ng , Ng and are the shape function and 0 ,S and w are the nodal variables. 

Using the notation of Eqn. (2.6) 3qn. (2.1) can be rewritten in the matrix form as 

L'^DL0 + S-^0 = O 
12 

-0-CiS + Vw = O (2.7) 

V^S + q-phw = 0 


where L, V are operators, Ci is matrix. 
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c “ 

' Oh 0 1 


dy dx 


and D is property matrix as 

, 1 V 0 

d=_1!i1^v 1 0 

12(l-v=) ^ ^ ^ 

L 2 J 

Now, using Galerkin method of deriving FEM equations with wei^ting functions 
Nq.Ns ,N^ on Equation 2.7 and integrating over the domain of the problem, Q, the 
following equations are obtained (Zienkiewcz, 1991, [22]) 

Ae + BS+E0 = f9 

B^0-PS + Cw = O (2.9) 


where 


CS + Fw = f^ 

A = j(LNe)DLNe^dn 

Q 

B = -jNeXdn 

a 

C = jN/(VN.)dn 

Q 

P = -jNs"C,Nsdn 

Q 

f.=jN,'Tdr 

r 


(2.10) 


E = ^jN,X<i« 
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F = phjN/N„dn 

o 

Here, the vector T represents the two moment components imposed on the boundary by the 
traction and F represents the imposed shear force. The shape functions Ne and Nw are chosen 
to have Co continuity but Ns can be discontinuous between elements. 

Indeed, such a choice allows S to be locally defined for each element and thus be 
e liminat ed at element level when P^K) and the main problem includes only the variables 0 
and w which have inertia components. 

Thus, equation (2.9) can be simplified by eliminating S (using the second equation) 
in terms of 0 and w, as; 


'a+BP''B^ BP“*cTe1 [E 0l 0 ^Tfe 

_ C^P'*B^ C^P''CjLwJ^Lo fJw 


( 2 . 10 ) 


in general this equation will be referred as , 

[K]'(Q} + [M]‘{Q} = (Fr (211) 

Such element equations are assembled to form the global equations as 

(I[k]’){Q} + (S[M]*){Q} = I{F}‘ 

The final set of assembled equations can be written as, 

where 

[m] = Yj [^]* > Global mass matrix 

, Global stiffcess matrix 
^ {f}* , Global applied force vector 

and {Ff^are Elemental Stif&iess Matrix, Elemental Mass Matrix and 

Elemental Force Vector. 


2.2.2 Integration Algorithm 

Eqtt (2. 12) represent a system of linear differential equation of second order and the solution 
to the equation can be obtained by standard procedures of solving differential equations. The 
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procedures for solving differential equations are divided into direct integration and mode 
superposition of which former are preferred in wave propagation probleuL In direct 
integration scheme, there are many different methods, which can be classified as ‘Explicit’ or 
‘Implicit’ whose advantages and disadvantages are given by (Bathe, 1990, [24]). 

In the present work, Newmark’s method of integration is used. 

Newmark’s method is an extension of ‘Linear acceleration method’, in which, linear 
variation of acceleration is assumed fi-om time t to / + A/ . 


Velocity and acceleration approximated in terms of displacements are 


t+At 


t+At 




— a 

j 




At"' 


(2.13) 

(2.14) 


where a and <5 ’ are the parameters chosen have accuracy and stability. 

Usually, a and 5 are taken as 0.25 and 0.5 respectively to get unconditional stability 

2.3 Energy Release Rate and Its Determination 

In fracture mechanics a crack can be characterized by four different parameters 

i. Energy release rate ( ) is energy based and is appUed to brittle and 

less ductile material 

ii. Stress intensity 6ctor ) is stress based This parameter is applied 

to brittle and less ductile material. 

iii. J-Integral (J), which has been developed to deal Avith ductile material and can 
also be applied to brittle material as well. 

iv. Crack tip opening displacement (CTOD) is displacement based and is 
developed for ductile materials. 

From computational point of view, stress intensity faetor requite special element so 
as to simulate stress singularity at the crack tip whereas J-Integral and energy release do not 
require any such modeling of stress singularity. In energy release rate approach, energy is 
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found out in entire domain and in J-Integral, evaluation is done along a path far away from 
crack tip. Thus both these approaches cleverly avoid analysis close to the crack tip and don’t 
need any modeling of stress singularity. 

Griffith (1920, [3]) outlined that in brittle fracture, as the crack extends new surfaces 
are formed. Formation of new surface requires surface energy. Consequently, crack in a 
brittle solid should advance when the reduction of the total potential energy of the body 
during a small amount of crack advance equals the surface energy of the new surface thereby 
created. Most of the energy release, as the crack advances, comes from the jmrts of the 
component, which are adjacent to the cracked surface. Two important parameters need to be 
considered 

1) How much energy is released when a crack advances 

2) Minimum energy required for crack to advance in forming two new surfaces. 

First parameter is measured in terms of energy release rate denoted by G. G is a 

function of crack size in general. Thus the energy release rate is defined as energy release per 
unit area extension of crack growth if there is a virtual crack growth, i.e. energy equal to G 
would be released from the system per unit extension of area. This energy is supplied by the 

elastic energy in the body and by the loading system. 

The energy requirement for a crack to grow per unit area extension is called as crack 

resistance and is usually denoted by symbol R. Crack resistance is sum of energy required to. 

1) Form new surface 

2) Cause anelastic deformation 

Both the available energy release rate and crack resistance are important to study the 
possibility of crack becoming critical. Obviously, when the available energy release rate far 
exceeds the crack resistance, the crack starts to grow at high speed. In the present study 
‘energy release rate’ is adopted, as it is a more comprehensive concept 

As the crack advances, 

1) Stifftiess of the component decreases. 

2) Strain energy in the component either increases or decreases. 

3) Work is done on the component if there is an application of load. 

4) Energy is being consumed to create two new surfaces. 
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Formulation for energy release rate is carried out invoking the principle of conservation 
of energy. First, consider thesis quasi-static crack growth case. Let the incremental increase in 
the crack area be AA. This crack growth is achieved by an incremental external work done by 
the external forces and the strain energy within the body of the component increases 

Let the incremental work by external force be AW^, change in strain energy be AU, 
and the available energy, GAA, satisfy energy balance equation as follows; 


GAA = AW^-AU 

(2.15) 

0--£(U-W„) 

(2.16) 

^ dn 


dA 

(2.17) 


where IT is the potential energy 

For a plate of uniform thickness B, dA= Bda 

where da is the incremental crack length, B is the thickness of the plate 


G = - 


1 dn 

B da 


(2.18) 


For dynamic crack propagation problem, the kinetic energy, T, of the body should be 
taken into consideration. Dynamic energy release rate, Gn, is different as some energy may 
be consumed to impart kinetic energy to the cracked piortion of the body and to generate 
stress waves. 

For dynamic case, energy balance becomes 

GM = AW^-AU-^ (2.19) 

where AT is the increment in kinetic energy in the body. 

For constant velocity crack propagation, 

G = ——(W^-U-T) (2.20) 

B da 


For crack moving with velocity v , 
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( 2 . 21 ) 



B 


In the present study, energy release rate is found out using equation (2. 19).Different energies 
are found out using 



(2.22) 

[uh\iQY[KM 

(2.23) 


(2.24) 

where, 


is the external work done 

U is the strain energy in the component 

T is the kinetic energy in the component 
{f} is the external force vector 
{(2}is the displacement vector 
^}is the velocity vector 
[at] is the global stififiiess matrix 
[M] is the global mass matrix 



2.4 Closure 

In this chapter, formulation of FEM equation is described along with integration 
algorithms for solving equilibrium equation governing linear dynamic response of a system 
of finite element, which will be used in the next chapters to construct a propagation element 
This chapter also describes some general concepts in energy and energy release rate. 
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CHAPTERS 


CRACK PROPAGATION ANALYSIS 


3.1 Introduction 

Over the last few years, physicists interested in nonlinear phenomena have returned to 
problems of crack propagation, and to the study of model problems in which cracks 
develop from the breaking of bonds. Also, it is now possible to model the entire 
development of a crack, in a fimte-element model of a continuum, even with inelastic 
effects (that is, plasticity) allowed for. Damage mechanics is well developed at the 
phenomenological level, and can also be built into codes that predict fracture. The 
admission of any model for microscopic processes introduces a length scale, and size 
effects may emerge. These are of practical importance as well as theoretical interest. 
Mathematical theory based on variational methods is developing, at least for cracks in 
elastic media. However, many problems remain; there is no adequate theory for the 
stability of a propagating crack, even in the frameworic of linear elasticity, at the 
present time. 

It is intended here to develop and improve fracture model (Stiffiiess Release 
Model, [19,20]) for simulating high-speed crack. 


3.2 Methods for Simulation of Dynamic Crack Propagation 
by Finite Element Method 

To simulate a crack propagation in solids two different concepts of finite element 
modeling are in use i.e., stationary mesh procedure and moving mesh procedure. Out 
of these two, moving mesh procedure, as the name implies involve change of mesh in 
each step as the crack propagates. This is computationally very intensive. Stationary 
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mesh procedures do not require changes in the mesh in general, and require only 
change in the boundary and loading conditions. The two procedures are explained in 
the following sections. 

3.2.1 Stationary Mesh Procedure 

In the simple stationary mesh procedure of mode ling Unear elastodynamic crack 
propagation, the nodes ahead of the crack tip are spaced at VcAt (Vc being crack 
velocity) and crack propagation is simulated by releasing the nodes one at a time. 
However, if a simple node release technique is used, release of constraint on the 
displacement (at the preceding crack tip) in each time step induces spurious high 
frequency oscillations in the finite element solution To overcome these dilBBculties, 
several algorithms have been suggested in the literature to release the node gradually 
over few time steps. 

Keegstra et al. [25, 26] si^ested a model in which node was not released 
instantaneously vdien the force at the crack tip reaches a critical value Fd, which is 
proportional to dynamic fracture toughness Kd. Instead the force was reduced from 
Fd gradually in accordance with an elastic spring model so that the nodal force 
vanishes when critical displacement is reached. Another node release techmque of 
gradual release of nodes is ‘Force release model’ vliich is described below. 

3.2.1.1 Force Release Model 

Suppose that the actual tip is located at ‘C’ in between the finite element nodes B and 
D as shown in Figure 3.1. The length segments BC and BD are b and d respectively. 
The holding back force, F, at node B is gradually reduced to zero over a number of 
time steps as the crack tip reaches to the node D. Various schemes available to 
decrease the force to zero are as follows; 

(i) Malluck and King [27] suggested the release rate based on constant stress 
intensity factor 
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where Fo is the original reaction force when the crack tip was located at node B. 

(ii) Rydhom et al. [28] Suggested the release rate based on constant energy 
release rate 




Kobayasi et al. 29] suggested the linear release rate based on no physical 
argument other that pure intuition. 


— = fl-- 


In order to have more gradual and smooth propagation of crack Kishore, Kumar 
Verma [30] used a modified method. Holding back force at the crack tip B is linearly 
decreased to zero when the crack reaches the end of the next element Thus when 
crack tip goes beyond node B then 


Fhb I d + d,} 

where is the reaction at node B, when the node was closed, b is the crack 
extension and d,di are the length of the elements as shown in Figure 3.1. And when 
the crack moves beyond node D to a point Di 




d + d, 


= 1 - 


{ d:+dj 

where Fhd is the reaction at node D, when the node was closed, b, b 2 are the crack 
extension and d, di, di are the element lengths as shown in Figure 3.1 
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In force release model, time increment At is usually taken such that it takes 15 
to 20 iterations for the crack tip to cross one element The amount of force necessary 
to keep the nodes together is determined and a factor is used to proportionately 
decrease the force as the crack tip advances. Thus, at subsequent times, though the 
problem is highly dynamic, the current calculations were based on the force 
calculated to keep the nodes together 15 to 20 time steps before. This will cause 
discrepancy as the crack advances further from one element to the next one causing 
large oscillations in the solution. This discrepancy should be overcome for obtaining 
better solutioiL A new model is proposed to overcome the drawbacks of force release 
model that instead of considering dynamic force at the crack tip takes into account the 
stiffiiess and mass of the element at the crack tip. It is described in detail in the later 
part of this chapter. 



Figure 3.1 Crack opening scheme in force release model 
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3.2.2 Moving Mesh Procedure 

In moving mesh procedure, entire mesh moves with the crack tip. It can be further 
classified into an approach wherein a) the entire mesh moves b) only the mesh in a 
finite and small region around the crack-tip moves along with the crack-tip. Aberson 
et al. [31, 32] used a singular crack-tip element as shown in Figure 3.2(a), which 
incorporates the first 13 terms of Williams’ eigenfunctions [33], appropriate to a 
stationary crack in a linear elastic body. The crack tip moves within the singular 
element, between the nodes A and B as shown in Figure 3.2(a). When the crack tip 
reaches the node B, the mesh pattern is changed suddenly, as illustrated in the figure. 



Figure 3.2 Discontinuously moving singular element 
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Another attempt at using Williams’ eigenfunctions was made by Patterson and 
Oldale [34, 35], The singular element (Figure 3.2(b)) has 13 nodes and is 
topologically equivalent to two assembled 8 noded isoparametric elements. The 
location of the sii^ular element is suddenly changed by a distance equal to the si 2 e of 
the regular element ahead of the crack tip, when the crack tip propagates a critical 
distance within the singular element. The displacements of the singular element are 
matched to those of the surrounding solid only at the nodes connected to the regular 
element. Thus, this model violates the displacement compatibility condition at the 
interfaces between the singular element and the surrounding regular elements. In 
another procedure, the singular element translates in each time step for which crack 
growth occurs. Thus crack tip is always at the center of the singular element 
throughout the analysis. The regular elements surrounding the moving singular 
element are continuously distorted. Mesh pattern around the singular element is 
periodically readjusted for simulating large amount of crack propagation. 


3.3 Stiffness Release Model for Crack Propagation 

The force release models to simulate crack propagation as applied to high-speed crack 
propagation has the following drawbacks. 

1 They give rise to very high oscillations in the solution. 

2. They are not suitable for solving inverse problem in crack propagation. If the 
crack propagates and stops in between an element and starts again later when 
another stress wave arrives, how much holding force to be used is not known. 

The time increment At, is usuaUy taken such that it takes 15 to 20 iterations 
for the crack-tip to cross one element Ihe amount of force necessary to keep the 
nodes together is determined and a factor is used to proportionaly decrease this force 
as the crack advances. Thus, at subsequent times, though the problem is highly 
dynamic, the current calculations were based on the force calculated to keep the 
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nodes closed 15 to 20 time steps before. This will cause discrepancy as the crack 
advances further results in large oscillations in the solution. 

In the case of decreasing energy release rate and possible crack arrest in between 
the nodes, using the force at particular value in the propagation problem will lead to a 
non-acceptable solutioiL This aspect acquires importance in the case of inverse 
problems. 



(c) 

Figure 3.3 Crack opening scheme in stiffness release model 
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In the present woit, a propagation model with modified stiffiiess and mass 
matrices is presented to take into accoimt the extension of crack in plates. In stiffness 
release model [19, 20], one-dimensional elastic element is added at the crack tip node, 
whose stiffness is reduced gradually as crack tip advances to the next node (Fig 3.3). 
Additional stiffiiess required is very large when crack tip is at the beginning of 
element and zero wben crack tip reaches to the other end. Amount of additional 
stiffiiess required is a function of crack length and is obtained fi'om quasi-static crack 
propagation analysis. 

Figure 3.3(a) shows a symmetric part of a deformed plate in which crack is up to 
the node B. Figure 3.3(c) shows the configuration of the plate vdien the crack has 
completely propagated to the next node C. Figure 3.3(b) shows the mocklling of the 
plate with partially propagated crack (some where between B and Q by the addition of a 
one-dimensional linear spring element whose stiffiiess is K^. 


be the stiffiiess at node B of original plate without additional spring 
Uq be the displacement of node when crack tip reaches the next node and 

element opens up completely. 

a: , can be found out firom equilibrium equation as 

Ko(u.-u) = K.u ( 3 - 7 ) 


K,^K, 



(3.8) 


Energy in the spring is given as 






(3.9) 


Energy release rate is written as 

^ AE (3.10) 

where AA is the increment in crack area 
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Let a non-dimensional parameter for crack length (a) be defined as, 
a=crack length (within element) / d 
where d is the distance between two nodes. 

AA = BdAa 

where B is the thickness of plate. 

For static case, energy release rate can be written as 


G = 



1 du 
Bd da 


Assuming linear variation of G and a, G can be written as 


G - C, +C2fl 


where, Q and Cj are the constants to be determined. 
Substituting G in Equation (3.11) 


^^1 


+ C2^ ~ ~ 


2 Bd 


f— 1 


(3.11) 

(3.12) 


(3.13) 


On integrating Equation (3.13), we get 


« = — 




Bd 


C^a + 


C,a' 


+ Cn 


(3.14) 


On substituting initial condition u=0 when a=0 and 
end condition u=uo vdien a= 1 


«o 


2 

*o“o 


Bd] 


C '\ 

^1+ — 

. 2 J 


Equation (3. 14) and Equation (3.15) gives 


(3.15) 
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(3.16) 


u _ 2C,a + 

Uq 2Cj +C 2 


On simplification it reduces to 


where 


— = a + C{a-a^) = f 
“0 



2C, +C 


2 2 


(3.17) 


(3.18) 


Equation (3.8) finally can be written as 




IzL] 

f J 


\Ndiere f=u/uo as given by Equation (3. 17) 


(3.19) 


Static stif&ess Kg can be determined for a given specimen under static condition. 

Energy release rate is determined using Equation (3.11) and (3.17). 

However as Equation 3.11 was based on quasi-static crack propagation, it 
doesn’t take into account the kinetic energy in the body. In dynamic problem total 
energy comprise strain energy in the body, kinetic energy and work done by external 
forces. For finding energy release rate, drop in the total energy in the system is 
considered in the present model and energy release rate is determined using 
fundamental Equation 2.19. 


3.4 Effect of Constant ‘C’ 

As can be seen, Ko has to be determined for the given plate under quasi-static 
conditions and given boundary conditions. However, it may be mferred that the near 
boundaries affect the value of Ko more than the far off boundaries. Suitable value of 
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'C can be determined by trial and error to obtain continuous energy release rate in 
between elements. For a quasi-static case usually it is a very small value. In the 
present problem, it is of the order of 0.25. 

However, it was observed that C has considerable effect on the dynamic crack 
propagatioiL Effect of change in value of C on f vs. non-dimensional crack length, a 
is shown in Figure (3.4) 

Spring stiflOness, K* that is gradually reduced as the crack propagates, is a 
function of crack length ‘a‘ and constant C as shown in Figure (3.5). Figure (3.5) 
shows that C plays major role in controlling the value of K®. Also constant, C is 
affected by change in element size and time step At For a given problem, its value 
needs to be chosen so that final energy release rate curve are smooth without wide 
fluctuatioBL Although the proposed stifhiess K* of the additional elastic spring 
element varies fi-om oc to 0, it changes discretely in a finite number of steps as the 
crack propagates through one element 
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f=a+C(a-a^) 


Ko* 


K«=K<,(l-f)/f 



Fignre 3.4 f vs. a curve for different C 



°01) 02 0.4 D^ flfi 10 

Non Dimensional crack length, a 


Figure 3.5 K, vs. a, for different C 
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3.5 Improved Stiffness Release Model 

Figure 3.5 shows that the variation in the stiffiiess, depends on the value of ‘C’, 
which also affects the energy release rate, G. The fluctuation in the energy release 
rate at element boundaries changes with the values of ‘C’. In order to have less 
fluctuation at the element boundaries one has to play with different values of ‘C’ and 
by trial and error method ‘C’ is chosen which has least fluctuatiorLTo overcome this 
difficulty, improved stif&iess release model is proposed as; 


f =a + C(a-a'')“ 


where a is an exponent 


K,=K, 


/ 


(3.20) 


(3.21) 


For a=l. Equation 3.20 represents original s tiffiie ss release model [19, 20]. Primary 
requirements for spring stiffiiess are as follows; 

Ks-> 00 as a -)• 0 and KgK) at a = 1 . The effect of the parameter a on the stiffiiess 
release model as well as on the energy release rate will be discussed in the next 
chapter. 


3.6 Mass Modification 


In the stif&iess release model [19, 20], as crack propagates, stiffiiess of the element 
decreases but mass matrix of the element remains constant In reality when the crack 
is at the beginning of the element, there is no mass to get accelerated and the effective 
mass increases gradually as the crack reaches to the end of the element At tiie end of 
the element, it becomes equal to the mass of the element To account for this shape 
functions of the elemental mass matrix at the crack tip are modified so that nodal 
masses vary as a fimction of non-dimensional crack length, a. 
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3.6.1 Modification of Shape Function and the Crack Tip Mass 
Matrix 

Consider the element at the crack tip as shown in Figure 3.6 with the crack partially 

propagated by an amount ‘a’ within the element. 

3 



Figure 3.6 Representation-of an element at the crack tip 


For a 9-noded triangular element the shape functions are given as: 


Ni =2L^i-L, 

N^=2I3,-L^ 


Nj =24-L3 


N,=4L,L3 


N,=4L3L, 

N 7 =-32L^, L 2 L 3 

Ng=-32L,L"3L3 
N, =-32L,L2L"3 


The variation of the shape function, is shown in Figure 3.7 


(3.22) 
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Figures.? Shape Function, Ni vs. Li 


The shape function, of the crack tip node is proposed to be a function as. 


N.=(2L^-L,/ 


(3.23) 


>^diere, 

pO(l.l-a) T 

a , a is non-dimensional crack length and 

r is an index. The function P is chosen such that P = 1 at a=l and very hi^ 
-> 00 as a -> 0 ] at beginning of the element as shown in Fig 3.8. 
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3D 



a 

Figure 3.8 p vs. Non-dimeusional Crack Length, a 


Thus, as the crack advances, the nature of the shape function, Ni varies in 
such a way that mass of the element increases gradually. Fig 3.9 shows variation 
of Ni at different values of ‘a’ which shows that area under Ni increases from 
zero (at a=0) to maximirm final value 1/6 (at a=l). The proposed shape function 
Ni along with p is justified at the end of this chapter. 
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Figure 3.9 Modified Shape, Ni vs. Li for difierent values of ‘a’ 


3.6.1.1 Modified Mass Matrix, [M]: 



0 

Nf 0 


0 N3 0 N4 

Nj 0 N 3 0.. 


N, 

0 


0 

N, 


(2x1 




[Rmi](i8xi8)=[N,f[NJ 


(3.26) 


Similarly, 

K]=k^ N, N, N, (3.27) 

[j^\6.6) = [N^ f IN^ ] (3.28) 


where, N 9 and Nw are shape functions for the nodal variables 0 and w 


respectively. 


CK*) 


where. 


[F\^,=ph\NjN,£l 

qC) 


p is density of the material, 
h is thickness of the plate. 


(3.29) 

(3.30) 


Finally modified mass matrix is given as [Section 2.2], 


[RMf ^ = 



0 


0 



(3.31) 


Equation 3.24 shows tiiat P depends on the index r. Figure 3.10 shows how P 
varies with index r and Figure 3.11 shows the variation of RM (1,1) as function of 
non-dimensional crack length, a, within the element From Figure 3.11 it can be seen 
that for r=0.9, the variation in element mass matrix is smooth. 
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By modifying shape function of elemental mass matrix, mass of the element is 
made to vary as a function of crack length such that the element mass is zero when 
the crack is at the beginning of the element and increases gradually till crack reaches 
to the end of the element, when mass of the element becomes equal to the mass of a 
regular element. Effect of this modification on energy release rate variation will be 
shown in chapter 4. 

3.7 Prediction of Crack Initiation, Propagation and Arrest 

The present method can be tqjplied to predict the crack propagation history under the 
known applied loads, if crack propagation resistance of the material is knowiL That 
means that fracture resistance of the material, G as a function of crack propagation 
velocity is supposed to be a known input data. At each time step, the available energy 
release rate is determined and it is compared with the material fincture resistance data 
to predict the crack propagation velocity. This is used to determine the crack 
propagation length in the given time interval and propagated accordingly. This 
procedure continues till crack arrest takes place. 


3.8 Closure 

In this chapter various methods of simulating dynamic crack propagation with 
stationary and moving mesh are described. The variation of various aspects of the 
proposed model is also discussed in detail. This model will be investigated with 
hypothetical data in ch^ 4. The effect of various parameters on energy release rate 
will be described in the succeeding chapter. 


40 



0 !o □ "2 nTi n~y« q *g 77a 

Figure 3.10 p vs. Non-dimensional crack length, a for different valnes of r 


Rine(l,l) [first diagonal ElenKut of Mass Matrix] vs. 
Iteration No for First Opening Element 



Iteration No 


Figure 3.11 RM (1,1) vs. Step No for different values of r 
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CHAPTER 4 

RESULTS AND DISCUSSION 


4.1 Introduction 

This chapter is concerned with the validation of the method and demonstrate the effect of 
improved stiffiiess release model along with mass modification on energy release rate and 
other fracture parameters. First FEM code developed is validated for static case. Later 
problems with a typical data are considered to discuss the effect of crack propagation 
model parameters along with mass modification on energy release rate, G. Finally a plate 
with a crack is analysed as: 

1 . Forward problem and 

2. Inverse problem 

in which plate subjected to impact loading is analysed by using the 9-noded hybrid 
triangular elements. 

4.2 Validation for Static case 

A plate without crack with all the four edges simply supported subjected to a uniformly 
distributed loading is analysed by the FEM code with the 9-noded hybrid triangular 
element and the results are compared with theoretical values. 

Dimensions of the plate; Length, L= 3 m and width, W= 2 m. 

The final refined mesh used is shown in Figure 4. 1 . 

The maviTmim deflection of the plate at its centre is given by the formula (Timoshenko, 
1959, [24]) 

Wmax=(tzqa'^)/D 

where,D = E*hVl2(l-.9^) 

E= Young’s Modulus of Elasticity 

h = thickness of the plate 
S = Poisson’s ratio 
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a = Numerical factor depending on the ratio b/a of the sides of the plate (for a 
uniformly loaded plate with the 4-edges simply supported) 

= 0.00772 

q = uniformly distributed load 
a = length of plate in x-direction= 2 m 
b = length of plate in y-diTection= 3 m 

For a steel plate, with E = 200 GPa, h = 0.03 m, 3 = 0.3, q = 1000 nW, a= 2 m , the 
above formula gives Wmax = 2.49 xlO"^ m The displacement contour obtained by FEM 
code is shown Figure 4.2. The value of maximum deflection Wmiiv obtained from the code 
is 2.6x 10"^ m. It can be observed that the error is 4.4 % for the 54-element mesh. 


4.3 Input Parameters for Hypothetical Problems 

In this section, the effects of various parameters in the crack propagation model are 
investigated for the plate subjected to typical load. The material properties, plate 
geometry, loading, and boundary conditions used for analysis are as follows. 

Material properties: 

• Modulus ofElasticity,E = 210 GPa 

• Poisson’s Ratio, v= 0.3 

• Density of the material, p = 7800 kg/m^ 

The rectangular plate is shown in Figure 4.3 has a central crack parallel to Y-axis. 
Dimension of the plate: 

• Length of the plate, L= 4 m 

• Width of the plate, W= 2 m 

• Crack Length, (2a )= 0.4 m 

• Thickness of the plate, h= 0.024 m 
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Figure 4.1 Meshing used for static case 



Figure 42 Displacement Contours for static case 
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Input Pulse: 

The plate is subjected to a uniformly distributed impulse load of intensity, q, which rises 
sharply and dies down slowly as shown in Figure 4.4 is given by, 

ampq__^ 

0.3679 

t =-^ 

*'fft max A 

4 

where tat_max= Rise time, 

pdur = Total pulse time 
Details of the input pulse are as follows: 

• Amplitude of the Load = 1000 N 

• Rise Time = 35 ps 

The plate has all the four sides simply supported is analysed as shown in Figure 4.5. As 
the problem is symmetric about X and Y-axis, only a quarter of the plate as shown in 
Figure 4.6, is analysed. 

Details of meshing as shown in Figure 4.7 is as follows: 

• No. of elements= 143 

• No. of nodes= 739 

• No. of degree of freedom= 1788 

Crack Propagation Details: 

• No. of element to open = 3 

• Time step (dt) = 2 ps 

• Load steps: Load steps are the number of steps to open one particular element. 
Load steps for first opening element = 400 

Load steps for second opening element = 240 
Load steps for third opening element = 220 
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Simply Supported edge, w=0 



Simply 
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Simply Supported edge, w=0 

Figure 4.5 Boundary Conditions 
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Figure 4. 6 Boundary Conditions for Symmetric Quarter plate 
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Figure 4.7 Mesh used for Hypothetical Problem 
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4.4 Results for Improved Stiffness Release Model 

4.4.1 Effect of the Parameter a on Spring Stiffness, Ks 

In the stiffness release model developed by Kishore, Reddy and Deshmukh [19, 20], the 
spring stiffness. Kg depends highly on the values of constant ‘C’ as seen from Figure 3.5. 
Figure 4.8 shows the variation in energy release rate for different values of ‘C’. ‘C’ plays 
major role in fluctuations of G at element boundaries. Figure 4.9 auH Figure 4.10 show 
fluctuations at the elemental boundaries. To overcome this difficulty, improved stiffiiess 
release model is proposed [Section 3.5] as 

f* = a + C(a-a^)“ and 


Figure 4. 1 1 to Figure 4. 16 shows the variation of Kg vs. non-dimensional crack length, a, 
for different values of a. For a=l nature of Kg largely depends on values of ‘C’. It is 
evident from Figure 4.11 that Kg becomes discontinuous at a= 0.5 for C= -2 as f 
approaches zero. Also for very high values of ‘C’ (C> 2), the stiffiiess of the spring is 
completely released \^en the crack is at the beginmng of the element itself, which is not 
realistic with the nature as stiffiiess. Kg should release over the whole element For 
example at a= 0.5 for C= 2, Kg becomes zero (Figure 4. 1 1 ). 

It is clearly observed that by increasing the value of ot, the dependence of spring 
stiffiiess. Kg on ‘C’ value reduces. From Figure 4.13, it can be seen that for a= 2, the 
range of ‘C’ value has been increased as compared to oc=l (Figure 4.11) for v4uch the 
propagation model works satisfactorily. Thus as the value of a increases, the range of ‘C’ 
for which model woiks satisfactorily increases. Hence with increase in oc, dependence of 
Kg on ‘C’ reduces. 
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Energy Release Rate, 




Figure 4.9 Energy Release Rate, G vs. No. of Iterations 


50 




Energy Release Rate, G 


Energy Rdtease Rate, G for DifTerent values of C for alpha= 1 (Interface betn 


element 1 and 2) 



Figure 4.10 Energy Release Rate, G vs. No. of Iterations 
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Figure 4.15 K, vs. a, for different C and ar=10 


Ko* 



Figure 4.16 K, vs. a, for different C and a=25 
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4.4.2 Effect of the parameter a on Energy Release Rate, G 

The effect of increase in a on energy release rate is shown in figure 4.17. It is observed 
that with increase in a, fluctuations in the energy release rate at the element boundaries 
decreases. The decrease in the fluctuations is fester initially for low values of a and 
becomes slower for higher values of a as seen fi'om the Figure 4.18. Thus moderate 
values of a between 2 to 10 are proposed for improved stiffiiess release model. 


4.5 Results for Mass Modification 

In stiffiiess release model [19, 20], ‘C’ value was chosen such that fluctuations in the 
energy release rate was Triinimum. In improved stiffiiess release model, the dependence of 
energy release rate on ‘C’ value as well as fluctuations at the elemental boundaries are 
reduced by introducing parameter a. Further by mass modification, the fluctuations at 

elemental boundaries can be stiU reduced. 

As the stiffiiess of elastic element decreases gradually to zero, the element mass 
can be increased fi-om zero to its usual value by modifying shape function Ni for an 
element containing propagating crack (Section 3.6.1). Energy release rate after this mass 
modification is shown in Figure 4. 19. Figure reveals that final energy release rate curve 
gives smooth variation within element and fairly stable solutions across element 

boundaries. 

Figure 4.20 shows energy release rate with and without mass modification. It can 
be seen very clearly that solution with mass modification is better t h a n one without mass 
modification (Figure 4.21 and Figure 4.22). Thus mass modification gives smooth 
variation in energy release rate within element as well as at the element boundaries. 
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Figure 4.17 Energy Release Rate, G vs. No. of Iterations 
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Figure 4.18 Energy Release Rate, G vs. No. of Iterations 




56 





bnergy Release Rate, G | Energy Release Rate, 



Figure 4.19 Energy Release Rate, G vs. No. of Iterations 
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Figure 4.20 Energy Release Rate, G vs. No. of Iterations 
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Figure 4.21 Energy Release Rate,G vs. No. of Iterations 


Energy Release rate, G with and without mass modification for 
alphsF 1 at the interface of the element 1 and 2 



—G without 


mass 

modification 


G with mass 
modification 


460 460 470 480 490 

No Iterations 


Figure 4.22 Energy Release Rate,G vs. No. of Iterations 
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Figure 423 Energy Release Rate, G vs. No. of Iterations 


Crack velocity vs. No. of iterations for Forward problem 
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Figure 4JJ7 Energy Release Rate vs. No. of Iterations 




Crack Length, a' 


Crack Length vs. No. of Kerations for Inverse problem 





CHAPTERS 


CONCLUSIONS AND SCOPE FOR FUTURE WORK 


5.1 Conclusion 


A modified stiflBtiess release model is proposed to simulate hi^-speed crack propagation. 
This model is further improved for mass modification and based on the results and 
discussion in Chapter 4, following conclusions can be drawn. 

1 . Parameter a reduces dependence of the sprii^ stiffiiess, Ks on C in stifhiess release 
model. 

2. For oo 1 , the fluctuations in energy release rate at element boundanes is reduced. 

3 . Modified model gives fairly smooth and more stable variation of energy release rate. 

4. Mass modification further reduces the fluctuations in energy release rate at the 
element boundaries 

5.2 Scope for Future Work 

1 . The crack propagation model may be fuither improved to get still better results. 

2. Present method can be extended to nonuniform crack speed. 

3 . This method can also be extended for 3-Dimensional crack propagation problems. 
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